1 



AN EXPLANATION FOR THE SOFT X-RAY EXCESS IN AGN 



J. Crummy \ A.C. Fabian', L. GaUo^, and R.R. Ross^ 

'institute of Astronomy, Madingley Road, Cambridge, UK. CB3 OHA 
^Max-Planck-Institut fiir extraterrestrische Physik, Postfach 1312, 85741 Garching, Germany 
^Physics Department, College of the Holy Cross, Worcester, MA 01610, USA 



ABSTRACT 



We present a large sample of type 1 AGN spectra taken 
with XMM-Newton. We fit them with the relativisti- 
cally blurred photoionized disc model of Ross & Fabian 
(2005). This model is based on an illuminated accretion 
disc of fluorescing and Compton-scattering gas, and in- 
cludes relativistic Doppler effects due to the rapid motion 
of the disc and general relativistic effects such as gravi- 
tational redshift due to presence of the black hole. The 
disc model successfully reproduces the X-ray continuum 
shape, including the soft excess, of all the sources. It pro- 
vides a natural explanation for the observation that the 
soft excess is at a constant temperature over a wide range 
of AGN properties. The model also reproduces many fea- 
tures that would conventionally be interpreted as absorp- 
tion edges. We use the model to measure properties of the 
AGN such as inclination, black hole rotation, and metal- 
Ucity. 
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clei; X-rays. 



1. INTRODUCTION 



The soft excess is an important component of the X-ray 
spectra of many AGN and is present in every source in 
this survey. The soft excess is defined as the enhanced 
emission below ^2 keV compared to an extrapolation of 
the approximately power law spectrum in the 2-10 keV 
band. This extra emission is generally approximately 
thermal in shape, and well fit with a black body of energy 
0.1 - 0.2 keV. This result stands over several decades 
in AGN mass, e.g. Walter & Fink (1993), Gierliriski & 
Done (2004), Porquet et al. (2004). This poses problems 
for the thermal interpretation of the data, as the tempera- 
ture does not scale in the expected way with luminosity 
(L oc T*), as well as being too high for the standard 
model of Shakura & Sunyaev (1973). Several alternative 
models have been proposed to account for one or both 
of these results, none of which have yet been broadly 



accepted by the community. This paper investigates 
photoionized emission blurred relativistically by motion 
in an accretion disc, which has been previously studied 
by e.g. Ballantyne, Iwasawa & Fabian (2001). We use 
the latest models from Ross & Fabian (2005), which 
include more ionization states and more recent atomic 
data than earlier versions (Ross & Fabian 1993). 



2. THE RELATIVISTICALLY BLURRED PHO- 
TOIONIZED DISC MODEL 



In the model of Ross & Fabian (2005) a semi-infinite slab 
of cold optically thick gas of constant density is illumi- 
nated by a power law, producing a Compton component 
and fluorescence lines from the ionized species in the gas. 
To produce the relativistically blurred photoionized disc 
model this reflected emission is added to the illuminating 
power law and the summed emission is convolved with 
a Laor (1991) profile to simulate the blurring from an 
accretion disc around a maximally rotating (Kerr) black 
hole. This assumes a geometrically thin and flat accretion 
disc with clearly defined inner and outer radii, and that 
the emissivity of the disc as a function of radius is 
described by a power law. This combined model allows 
fitting to the inclination of the disc, the inner and outer 
radii (the outer radius is generally poorly constrained as 
emission is centrally concentrated), the emissivity index 
of the disc, the iron abundance of the gas, the spectral 
index of the illuminating power law and the ionization 
parameter (which we define as ^ = AttF/uh, where F 
is the illuminating energy flux and uh is the hydrogen 
number density in the illuminated layer, a measure of the 
ratio of the energy density in the illuminating radiation to 
the atomic number density in the gas). Since the model 
is physically motivated measuring these free parameters 
can give us information about the sources. An example 
of the model is shown in Figuren(bottom right panel). 
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Figure 1. This figure shows a comparison between the thermal model and the relativistically blurred photoionized disc 
reflection model for NGC 4051. The top panels show spectra and residuals to the thermal (left) and disc reflection ( right) 
models, and the bottom panels show the components of the thermal (left) and disc reflection (right) models. Both include 
a power law (dotted line) and narrow iron emission (sharp line at ^6.4 keV), the components shown in dashed lines are 
a black body and the relativistically blurred disc reflection respectively. As the residuals show, the disc reflection model 
is a much better flt. Most of the remaining residuals to NGC 4051 in the soft band (top right panel, 0.5 - 1.0 keV) have 
been shown to be due to narrow line emission (Ponti et al., in prep). 



3. DATA 



We used publicly available archival XMM-Newton 
data on 22 type 1 AGN from the Palomar-Green (PG) 
sample and a selection of 12 other Seyfert 1 galaxies 
with high-quality observations available. The sample 
includes several well-studied AGN, such as NGC 4051, 
PG 1211+143 and I Zw 001. We analyse the longest 
available observation where the EPIC pn camera took 
data. We reduced the Observation Data Files in the 
standard way using SAS 6 . to produce spectra (See 
Crummy et al. 2005b for details), taking the range 0.3 
- 12.0 keV as the region in which the pn is accurately 
calibrated. We do not include MOS data in our fits. We 
grouped the spectra so each bin includes at least 20 
source counts, so statistics are applicable. Quoted 
errors are 90 per cent limits on one parameter (Ax^ = 
2.706). 



4. ANALYSIS 



We fit the data with two classes of model using xspec, 
the standard power law with a black body to model the 
soft excess, and the relativistically blurred photoionized 
disc reflection model. Both models are subject to 



absorption from cold gas in our Galaxy, we fix the 
amount of absorption at the value given by the nh ftool 
(Dickey & Lockman 1990). We also aUow for absorbing 
matter at the AGN by fitting for extra cold absorption 
and up to two absorption edges in the 0.45 - 1.1 keV 
band. These components are redshifted such that they 
are local to the AGN. We finally allow for cold reflection 
(e.g. from distant gas such as a torus) by including a 
narrow iron line at 6.4 keV (in the source frame). We 
are concerned primarily with the underlying continuum 
rather than these features, so we include in our final 
model only those components which improve the overall 
by > 2.7 per lost degree of freedom. Since we are 
performing fits over the entire XMM-Newton band, we 
only report lines for which xspec reports a non-zero 
minimum equivalent width, this avoids fitting spurious 
lines due to a curvature in the continuum not addressed 
by the model. This conservative procedure means faint 
lines are likely to be missed, although any line that 
significantly affects the overall goodness of fit will be 
included. In xspec terminology our models are phabs 

* zphabs * zedge * zedge * (powerlaw 

+ zbbody + zgauss and phabs * zphabs 

* zedge * zedge * (kdblur (powerlaw + 
atable{ref lion}) + zgauss), where kdblur 
is a convolution model using a Laor line and re f lion 
is a table model of reflection from cold gas (including a 
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Figure 2. Example disc reflection model flts plotted in vF^, showing the amount of flux emitted as a function of energy. 
The crosses are data, the data-following (blue) line is the fit, the smooth (green) line the power law component and the 
spiky ( red) line the reflection component. The inserts show the shape of the Laor line used in the convolution. The soft 
excess is shown to be the result of blurred line emission, and the extreme blurring effect on the iron lines around 6.4 keV 
is visible. The flgure also illustrates the effect of ionization parameter on reflection spectra, for PG 1309+355, = 3, PG 
1501+108, ^ = 510, PG 1116+215, ^ = 1270 and ARK 564, ^ ^ 3120. A colour version is available online. 



redshift)^ 



5. RESULTS & DISCUSSION 



Our results are presented in full in Crummy et al. 2005b. 
One of the sources in the sample, PG 1404+226, has 
been investigated with the same model in a previous pa- 
per (Crummy et al. 2005a). We find that in general the 
relativistically blurred photoionized disc reflection model 
is a better fit to the data, with 25 of the 34 sources show- 
ing an improvement in > 2.7 per degree of freedom 
(note that this does not correspond to 90 per cent proba- 
bility, the distribution is not calibrated across models). 
6 sources show a significant worsening and 3 sources are 
inconclusive. In some sources the improvement is very 
marked, see Figure^] 

The disc reflection model reproduces the shape of the 
continuum well, and naturally explains the constant tem- 
perature of the soft excess; see Figure|2] The model also 

'Models available from 
lhttp://www-xray.ast. cam. ac.uk/~jc/kdblur.html and 
|http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/models/reflion.html| 



reproduces many features that would be otherwise inter- 
preted as absorption edges in medium resolution pn data. 
When the thermal model is used, 17 of the 34 sources 
appear to have an absorption edge. This reduces to just 
7 sources when the disc reflection model is used. The 
disc reflection model is somewhat less smooth than the 
thermal model (see Figure and bumps on the scale 
of absorption edges are possible. Absorption lines can- 
not be reproduced by the model, therefore the possibil- 
ity of outflows from absorption edges should be checked 
by investigating reflection models, or by detecting ab- 
sorption lines from outflowing matter This has impli- 
cations for the energetics and cosmological evolution of 
these systems, as some measured outflows can involve 
similar amounts of power as the radiation e.g. Pounds et 
al. (2003). The model also explains the apparent lack of 
broad iron lines in most AGN spectra. Very few broad 
iron lines have been clearly detected (the most well stud- 
ied being MCG -6-30-15, Tanaka et al. 1995), which 
is unexpected since it is believed that these systems are 
powered by accreting matter near the black hole, where 
broad iron lines are formed. Our results suggest that the 
broad iron lines do occur, but are broadened to the point 
of near undetectability. A previous paper by Gallo et 
al. (2004) supports this hypothesis, they fit an extremely 
broad excess of counts in the hard band of MRK 0586 
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Figure 3. This histogram shows the measured inclina- 
tions of the sources (the solid line) and the expected 
distribution if the inclinations were random ( the dashed 
line). There is a deficit at high inclinations. 



with a Laor profile. MRK 0586 is included in our sample, 
and we find an inclination consistent with theirs; we even 
fit consistent inclinations when we exclude the hard band 
from the fit and determine the inclination purely from the 
blurring effects on the soft excess. This illustrates how 
the hard and soft bands are connected, they key to the ab- 
sence of obvious broad lines is in the shape of the soft 
excess. Conversely, it is clear that the soft excess is noth- 
ing more than a blend of many broad lines. 
As noted in the previous paragraph, it is also possible to 
measure the inclination of the central disc using the disc 
reflection model. Our measured inclinations are plotted 
in Figure|3] We find that the inclinations are inconsistent 
with being random, with a Kolmogorov-Smirnov (K-S) 
test probability of 0.06. Performing a K-S test against a 
distribution which is random over 0° - 81° (ignoring any 
data points above 81°) gives a probability of 0.34, so the 
data is somewhat consistent with being random over this 
range. The deficit at low inclinations may be a selection 
effect, discs at ~0° Doppler beam all of their emission 
along the plane of the disc and so would appear dim. The 
model also takes no account of limb darkening, which is 
strong at low inclinations. At high inclinations the deficit 
could be due to torus obscuration, with the sources that 
do have high inclination perhaps not having a torus, or 
with the torus unaligned with the central disc. 
Abundances may also be measured with the model. The 
model assumes a fixed, solar abundance for all elements 
except iron, for which the abundance is a free parame- 
ter. The measured abundance is consistent with being be- 
tween solar and 1/3 solar for 28 of the 34 sources. No 
sources have measured abundances below 1/3 solar, and 
6 have abundances which are clearly above solar. These 
fairly low iron abundances are another factor in the appar- 
ent absence of broad iron lines. The inability to fit more 
than one element abundance means these results are ten- 
tative, but show that in principle the composition of AGN 



accretion discs can be measured from their X-ray spectra. 
Finally, the rotation of the central black hole may be mea- 
sured using this model. The model uses convolution with 
a Laor line, which is based on a maximally rotating (Kerr) 
black hole. It is equally possible to use a model based on 
a non-rotating (Schwarzschild) black hole (diskline 
in xspec, Fabian et al. 1989). We also performed fits to 
all our sources using this non-rotating model, and found 
that in all cases the non-rotating model has a worse good- 
ness of fit, with only two of the 34 sources where the fit 
is of comparable quality. This shows that rotating black 
holes dominate our sample. The black hole rotation is 
not a free parameter in our analysis, but the inner radius 
of the accretion disc is. The last stable orbit around a 
non-rotating black hole is at 6 gravitational radii, and the 
last stable orbit around a maximally rotating black hole 
is at 1.235 gravitational radii. Inside this radius matter 
plunges into the black hole, and does not emit strongly 
(Section 3.4 of Fabian & Miniutti, 2005). All 34 sources 
have inner disc radii consistent with being below 6 gravi- 
tational radii, and 29 of the 34 have inner radii consistent 
with being below 1.3 gravitational radii. This is a strong 
indication that most black holes are maximally rotating, 
in agreement with theoretical predictions, e.g. Volonteri 
et al. (2005) predicts that 70 per cent of AGN black holes 
are maximally rotating. 

One further interesting result from the fits is that the illu- 
minating power-law component is often undetected. This 
seems to be a problem, as it is hard to imagine how the 
disc could be illuminated by radiation we cannot see. 
However, it is quite possible for relativistic effects due 
to the black hole to "hide" the illuminating continuum 
from us. If the source of the illuminating continuum is 
above the disc, the light from it will be bent and impact on 
the disc, with little escaping to be observed (Miniutti & 
Fabian 2004). Emission from the disc is Doppler beamed 
away from the black hole and will still escape. This pri- 
marily applies to sources which do not share the rotation 
velocity of the disc, e.g. the base of a weak jet or shocks 
in a failed jet (Ghiselhni, Haardt & Matt 2004). 



6. CONCLUSIONS 



We fit a large number of type 1 AGN with a thermal 
model and a relativistically blurred photoionized disc re- 
flection model based on Ross & Fabian (2005). We find 
that: 

• The disc reflection model fits the data better. 

• The disc reflection model reproduces all the major 
features of all the sources, including the soft ex- 
cess and the absence of obvious broad iron lines. 
The model explains the constant temperature of the 
soft excess across all the sources, since it is formed 
from a large number of highly broadened lines. The 
model reproduces many features that might other- 
wise be interpreted as absorption edges. 

• Black holes in AGN strongly rotate. 
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• The central discs around type 1 AGN have a wide 
distribution in inchnation, with possible deficits at 
very low and high incUnations. The high inclination 
deficit may indicate torus obscuration. 

• The elemental abundances in accretion discs can be 
measured, iron abundances in these sources are solar 
or mildly sub- solar. 

The relativistically blurred photoionized disc reflection 
model is an important tool in the study of AGN. Taking 
account of their intrinsically relativistic nature answers 
several questions about their spectra, as well as providing 
information about the central regions. 
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